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Abstract — In this paper, optimal linear transceiver designs 
for multi-hop amplify-and-forward (AF) Multiple-input Multiple- 
out (MIMO) relaying systems with Gaussian distributed chan- 
nel estimation errors are investigated. Some commonly used 
transceiver design criteria are unified into a single matrix-variate 
optimization problem. With novel applications of majorization 
theory and properties of matrix-variate function, the optimal 
structure of robust transceiver is first derived. Based on the 
optimal structure, the original transceiver design problems are 
reduced to much simpler problems with only scalar variables 
whose solutions are readily obtained by iterative water-filling 
algorithms. The performance advantages of the proposed robust 
designs are demonstrated by the simulation results. 

I. Introduction 

In order to satisfy the emerging requirements for high 
speed ubiquitous wireless communications, MIMO coopera- 
tive communication has become one of the key parts in the 
future wireless standards such as LTE, IMT-Advanced, Win- 
ner project, etc. Transceiver design for amplify-and-forward 
(AF) MIMO relaying systems has been reported in ITI-H. 
There are various design criteria with different goals. The 
most common criteria are capacity maximization H], 121, lEI 
and data mean-square-error (MSE) minimization [Sl-llSl. In 
most of the previous works on transceiver design, most of 
the designs are restricted for dual-hop relaying systems and 
furthermore channel state information (CSI) is assumed to be 
perfectly known. Unfortunately, channel estimation errors are 
inevitable in practical systems. To mitigate the effect on the 
performance of AF relaying systems, such channel estimation 
errors should be taken into account in the transceiver design 
process. 

In this paper, we consider robust transceiver design for a 
multi-hop AF relaying system with channel estimation errors. 
Taking the Gaussian distributed channel errors into account, 
the precoder at source, multiple forwarding matrices at all the 
relays and equalizer at destination are jointly designed. The 
structure of the optimal solution for the unified problem is 
derived based on Majorization theory and properties of vector- 
monotone functions. The derived optimal structure covers most 
of the existing transceiver design results in point-to-point 
and dual-hop AF MIMO relaying systems as special cases. 
With the optimal structure, iterative water-filling solutions 
are proposed to obtain the remaining unknown parameters in 



the transceiver The performance advantages of the proposed 
robust designs are demonstrated by simulation results. 

The following notations are used throughout this paper 
Boldface lowercase letters denote vectors, while boldface 
uppercase letters denote matrices. The notation denotes the 
Hermitian of the matrix Z, and Tr(Z) is the trace of the matrix 
Z. The symbol Im denotes the M x M identity matrix, while 
Om.n denotes the M x N all zero matrix. The notation Z^/^ is 
the Hermitian square root of the positive semidefinite matrix 
Z, such that Z^/^Z^/^ = Z and Z^/^ is also a Hermitian 
matrix. The symbol Ai(Z) represents the i**^ largest eigenvalue 
of Z. The symbol (g) denotes the Kronecker product. For two 
Hermitian matrices, C ^ D means that C — D is a positive 
semi-definite matrix. For two vectors, x > y represents that 
each element of x is larger than the corresponding counterpart 
of y. The symbol A \ represents a rectangular diagonal 
matrix with decreasing diagonal elements. 

II. System Model 

In this paper, a multi-hop AF MIMO relaying system is 
considered. There is one source with A^i antennas wants to 
communicate with the destination with Mk antennas through 
K — 1 relays. For the fc'^ relay, it has Mk receive antennas 
and Nk+i transmit antennas. It is obvious that the dual-hop 
AF MIMO relaying systems is one of its special cases when 
K = 2. 

At the source, a x 1 data vector s with covariance matrix 
Rs = E{ss^} = I^r is transmitted through a precoder matrix 
Pi. The received signal xi at the first relay is 

xi=HiPiS + ni (1) 

where Hi is the MIMO channel matrix between the source 
and the first relay, and rii is the additive Gaussian noise vector 
at the first relay with zero mean and covariance matrix R,ij = 

At the first relay, the received signal xi is first multiplied 
by a forwarding matrix P2 and then the resultant signal is 
transmitted to the second relay. The received signal X2 at the 
second relay is given by 

X2 = H2P2X1 + n2 

= H2P2HiPiS + H2P2ni+n2, (2) 



where H2 is the MIMO channel matrix between the first relay 
and the second relay, and n2 is the additive Gaussian noise 
vector at the second relay with zero mean and covariance 
matrix R„2 = crfljv/^. Similarly, the received signal at k*"^ 
relay can be written as 



Xfc = HfcPfeXfc_i +nfc 



(3) 



where Hfc is the channel for the k hop, and is the additive 
Gaussian noise with zero mean and covariance matrix R„ — 



Finally, for a K-hop AF MIMO relaying system, the re- 
ceived signal at the destination is 



K 



Lk=l 
K 




nK, (4) 



where nfc=i^fe denotes Z/^ x • • • x Zi. In order to guarantee 
the transmitted data s can be recovered at the destination, it is 
assumed that Nk and Alk are greater than or equal to ||3]. 

In practical systems, because of limited length of training 
sequences, channel estimation errors are inevitable. With chan- 
nel estimation errors, we can write 

Hfe = Hfc + AHfc, (5) 

where H;; is the estimated channel in the k*"^ hop and AH^ 
is the corresponding channel estimation error whose elements 
are zero mean Gaussian random variables. Moreover, the 
Mk X Nk matrix AH^ can be decomposed using the widely 
used Kronecker model AHfc = 'El^^Uw.k'^l^'^ Q, H). The 
elements of the Mk x Nk matrix llw,k are independent and 
identically distributed (i.i.d.) Gaussian random variables with 
zero mean and unit variance. The specific formulas of the 
row correlation matrix Sj, and the column correlation matrix 
S&fc are determined by the training sequences and channel 
estimators being used IQ, ID. 

At the destination, a linear equalizer G is employed to 
detect the desired data vector s. The resulting data MSE 
matrix equals to $(G) = E{(Gy — s)(Gy — s)^}, where 
the expectation is taken with respect to random data, channel 
estimation errors, and noise. Following a similar derivation in 
dual-hop systems |8J, the MSE matrix is derived to be 



*(G) 
= E{(Gy 



s)(Gy-s)H} 





" K 


H |- 


R-?iA-]G^ + Iw — 


n SfcPfc 


G" G 




.fe=i 





Tr(PKRxK-iPi*if)Sif 

K 

n SfePfc 



where the received signal covariance matrix Rx^, 
relay satisfies the following recursive formula 

Tr(PfcRx,_,pH*,)£ 



(6) 

at the fc*^ 



R 



HfcPfcRx, P^H" 



and Rxo 
at the source 



Rs 



(7) 

In represents the signal covariance matrix 



III. Transceiver Design Problems 

A. Objective Functions 

There are various performance metrics for transceiver de- 
signs. In the following, we focus on two widely used metrics. 

(1) In general, for balancing the performance across dif- 
ferent data streams, (e.g., minimizing the worst data stream 
MSE), the objective function is written as |i9j 

Obj 1: i,,[d{^{G))] (8) 

where ipii') is an increasing Schur-convex functiorQ and 
d(*(G)) = [[*(G)]ia [*(G)]2.2 ---f, with the symbol 
represents the entry of Z. 

(2) On the other hand, if a preference is given over a certain 
data streams, (e.g., loading more resources to the data streams 
with better channel state information), the objective function 
can be written as 

Obj 2: V'2[d(*(G))] (9) 
where ip2i') is an increasing Schur-concave function. 

B. Problem Formulation 

Although the above two criteria aim at different designs, the 
transceiver design optimization problem can be unified into a 
single form: 



min /(*(G)) 
s.t. Tr(PfeRx,_,Pf ) < Pk, 



k = 1, 



,K (10) 

where the objective function /(•) is a real-valued matrix- 
variate function with $(G) as its argument. Notice that for all 
the two objectives described above, /(•) is a matrix-monotone 
increasing function. 

For ([Tol l, there is no constraint on the equalizer G. We can 
differentiate the trace of (|6]l with respect to G and obtain the 
LMMSE equalizer 

H 

na.p. 

.k=l 

+ Tr(PKRxA-iP|*K)SA' + R„k] 
with the property [lOl 

*(Glmmse) ^ *(G). 

Because /(•) is a matrix-monotone increasing function, (fT2] | 
implies that Glmmse minimizes the objective function in dTol l. 
Substituting the optimal equalizer of ( fTTT i into $(G) in (|6]l, 
^•(G) equals to 

-I H 



'LMMSE 



[H/fPifRxK-iPi-H^ 



(11) 



(12) 



MSE — i-N 



K 



n HfePA 



.k=l 



[Ha'Pr-Rx 



pHfVH 



Tr(PxRxA-_iP|*K)S/^ + R„,,]- 



K 



(13) 



.k=l 



'The detailed introduction of Schur-concave/convext functions, and ma- 
jorization theory is given in II II . 



For multi-hop AF MIMO relaying systems, the received 
signal at k*"^ relay depends on the forwarding matrices at all 
preceding relays, making the power allocations at different 
relays couples with each other (as seen in the constrains of 
([Toll), and thus the problem ( fTOl i difficult to solve. In order to 
simplify the problem, we define the following new variable in 
terms of Pfc: 



Property 1: At the optimal value of P 1, must have the 
structure of 

= Uj,diag[A(0)]Uf^ (18) 

where the vector A(0) = [Ai(0), • ■ • , AAr(0)]T with A„(0) 
being the n*'' largest eigenvalue of 0, and 



— 1/2 



(14) 



njM^ and Qfc is an 



where Kp, = Tr(FfeFH*fe)Sfc + a 
unknown unitary matrix. The introduction of Qfc is due to that 
fact that for a positive semi-definite matrix M, its square roots 
has the form M^/^Q where Q is an unitary matrix. Notice 
that Fi = Pi. With the new variable, the MSE matrix $mse 
is reformulated as 



MSE 



K 



l/2j^-l/2 



HfeFfc 
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Af---AlAK---A,. 



(15) 



Meanwhile, with the new variables Ffc, the corresponding 
power constraint in the fc"^ hop can now be rewritten as 



Tr(FfcFH) < Pfc. 



(16) 



It is obvious that with the new variables Ffc, the constraints 
become independent of each other. Putting (fTsT i and (fT6] l into 
([Tol l, the transceiver design problem can be reformulated as 



P 1: min 
s.t. 



0) 



f{lN 

Tr(FfcFH)<Pfc, k 
& = Af---AlAK^ 

Qfc Qfc = iMfc 



= I,-- 

■A, 



,K 



(17) 



From the definition of Ak in (flST l and noticing that Kpj. = 
Tr(FfcFf *fc)£fc + aljM,, it can be seen that Ffc appears 
at multiple positions in the objective function. Therefore, the 
optimization problem is much more complicated than the 
counterpart with prefect CSI. Indeed, as demonstrated by 
existing works, robust transceiver design for point-to-point or 
dual-hop relaying MIMO systems is much more complicated 
and challenging than its counterpart with perfect CSI Q, (HJ. 

IV. Optimal Structure of Robust Transceiver 

Based on the formulations of the objectives given in (|8]l and 
(|9]l, in Appendix |A] it is proved that P 1 has the following 
property. 



Ur 



(19) 



Qf for Obj 1 
In for Obj 2 

In ( fT9l l. the unitary matrix Uw is defined from the eigen- 
decomposition W — Uw^wU^ with Aw \, the matrix 
UAib is an arbitrary unitary matrix, and Qp is the unitary 
matrix which makes QFdiag[A(0)]Qp having identical di- 
agonal elements. Furthermore, with this optimal structure, the 
objective function of P 1 equals to 

f{lN - Undiag[A(0)]Uf^) = 5f[A(0)] (20) 
& 

where g(») is a monotobically decreasing and Schur-concave 
function with respective to A(0). 
Proof: See Appendix |A] ■ 

Based on Property 1, the objective function of ( fTTI i can 
be directly replaced by g[A(0)] and thus the optimization 
problem is simplified as 



P 2: min g[A(0)] 

-T k ,'-ik 



S.t. = Af 

Tr(Ffc.FH) < Pfc, Q^Qfc-lA/, 
= Uodiag[A(0)]UH 



(21) 



where Afc's are defined in ( fTSl ). In order to further simplify 
the optimization problem, we make use of the following two 
additional properties. 

Property 2: As is a decreasing and Schur-concave 

function and A(0) 7(0), the objective function in P 
2 satisfies 



g(A(0))>g([7i(0) ■■■jNi®)!') 



(22) 



=-f{&) 



with 7,(0) ^ A,(A|Ak)A,(AH_iAk_i) • • • KiAfAi), 

(23) 

holds when the neighboring Afc's 



where the equality in 
satisfy 



Va, =Ua,_i, k 



,K 



(24) 



with unitary matrices Va^ ^"d V^j, being defined based 
on the following singular value decomposition Ak ~ 
Ua.Aa.VH^ with Aa, \. 

Property 3: As g{») is a monotonically decreasing function 
with respective to its vector argument, the optimal solutions 
of the optimization problem always occur on the boundary: 



Tr(FfcFH) = Pfc 



(25) 



-The specific expressions of g(») are given in Appendix lA] but they are 
not important for the derivation of the optimal structures. 



Furthermore, defining 



2 



(26) 



with ttfc = Tr(Sfc)/Mfc which is a constant, (l25i l is equivalent 
to 

Tr[FfeFH(afcP,*,. + <I^J]/,,^^. = P,. (27) 

Based on Properties 2 and 3 , the optimal solution of 
the optimization problem (l2Tl i is exactly the optimal solution 
of the following new optimization problem with different 
constraints 

P 3: min gh{®)] 

s.t. TT[FkF^{akPk^k + 'yllN,)]/rif, = Pk 
= 

Q^Qfc = Im, , = Uj,diag[7(0)]UH 
Va, =Ua,_,, k^2,---,K. (28) 

Noticing that gr(») is a monotonically decreasing function, 
solving P 3 gives the following structure for the optimal 
solution. 

Conclusion 1: Defining unitary matrices \J-Hk ^nd V-^^ 
based on the following singular value decomposition 

(KF,/7y/,)-i/2H,KP,.*,. + <72Ja.J-i/2 = U«^A„,VH^. 
with A-H,\ and U«„=Un, (29) 

when 'S'fe cx I or oc I, the optimal solutions of the 
optimization problem ( |28] | have the following structure 



i = 1, ■ • ■ , A^, the optimization problem for computing Ajr^, 
becomes 

min g[7(0)] 

N 

s.t. Y.fh-P^ 

1=1 

7(0) = [71(e) •••7iv(0)F 

''kA 



7.(e) 



\.\.k=lJk.i"'k 



nti(/i>L + i) 



(32) 



The specific methods for finding ^ depend on the ex- 
pressions of g(»). In the following, we discuss the solution 
of (I32I 1 in more detail. The design criterion of MAX-MSE 
minimization is taken as example to show how to compute 

MAX-MSE minimization is a special case of Obj 1 in ([Sll 
and in this case, ■0]^ (d($MSE)) = max[$MSE]i,i- Further- 
more, in Appendix lAl it is proved that g(A(0)) = ■0j[ljv — 
{T.^=iH®)/N) ® In]- Therefore, ff[7(0)] equals to 

g[7(0)] = max [1n - (^^^ j«(0)/iV) ® In 
1 ^ 



Ffc,opt = v/^,(A^J(afcPfcVi>, + ^2 

Qfc.opt = , (30) 

where Vfi^, jv and U>ij, jv are the matrices consisting of the 
first N columns of V^j. and U-^j., respectively, and A^^. is 
a N X N unknown diagonal matrix. The scalar ^^.(A^j.) is a 
function of A^F^and equals to 

^fc(A:Fj = Vf, 

= </{! - a^TiiV^^M^kPk^k + tIInJ-'/' 
X *,(a,F,*, + <jllNj-'/^Vn,,NA%J}. 

(31) 

In the optimal structure given by ( l30b . the scalar variable 
^j.(A:Ffc) is only a function of the matrix Aj^^ and therefore 
the only unknown variable in dSOl l is Aj^^. The remaining 
unknown diagonal elements of A^^r^^ can be obtained by water- 
filUng alike solution as discussed in the next section. 

V. Computations OF A IT,, 



based on which the optimization problem ( [32] ) becomes 

, 1 'n.k=lifk,i^'k.i) 

mm 1 > — 

N 



(34) 



The optimization problem (|34] i can be solved by using iterative 
water-filling algorithm. 

VI. Simulation Results and Discussions 

In this section, the performance of the proposed robust 
designs are evaluated by simulations. For the purpose of com- 
parison, the algorithms based on the estimated channel only 
(without taking the channel estimation errors into account) 
are also simulated. In the following, we consider a three- 
hop AF MIMO relaying system where all nodes are equipped 
with 4 antennas. Furthermore, the estimation error correlation 
matrices are chosen as the popular exponential model JU i.e., 
[^k]i.j = a'*^-'' and [Sfc],;.j = The estimated channels 

fife's, are generated based on the following complex Gaussian 
distributions 



Hi 



CM 



Ah,Nk{OAh.Nk 



(35) 



The remaining unknown variables in ( I30b are only Ayr^. 
Substituting the optimal structures given by Conclusion 1 into 



P 3 and defining [A^^Ji 



and [Ayrji^i = /fe^, for 



such that channel realizations Hfe = fife + AHfe have unit 
variance. We define the signal-to-noise ratio (SNR) for the 
k*"^ link as Pk/a'^^. At the source node, four independent 
data streams are transmitted and in each data stream, N^ata = 



♦ Capacity maximization witti estimated CSI only 
-V- Capacity maximization robust design 




SNR (dB) 



Fig. 1. BERs of the proposed robust design with different design objectives, 
when a = 0.4, ^ = and ct^ = 0.004. 

10000 independent QPSK symbols are transmitted. Each point 
in the following figure is an average of 10000 trials. 

Fig. [1] shows the bit error rate (BER) of the proposed robust 
designs with different performance metrics: sum MSE mini- 
mization, mutual information maximization and MAX-MSE 
minimization. Other parameters are taken as a = 0.4, /3 = 
and CTg = 0.004. It can be seen that the former two criteria 
have better performance than the latter one. Meanwhile, the 
capacity maximization based on estimated CSI only is given 
to show that the proposed robust designs are better than that 
of the design with estimated CSI only. 

VII. Conclusions 

Bayesian robust transceiver design for multi-hop AF MIMO 
relaying systems under channel estimation errors was consid- 
ered. Various transceiver design criteria were unified into a 
single optimization framework. Using majorization theory and 
properties of matrix-variate functions, the optimal structure of 
transceivers was derived. Then, the transceiver design prob- 
lems were greatly simphfied and the remaining unknowns were 
obtained by iterative water-filling solutions. The performance 
of the proposed transceiver designs has been demonstrated by 
simulation results. 

Appendix A 
Proof of Property 1 

Obj 1: For the diagonal elements of the positive semi-definite 
matrix $mse = Iat — 0, we have the following relationship 

mi 

EN 
K{@)IN)®1n <d{lN -@) (36) 
■i—i 

with the equality holds if and only if [0]i,i ~ J2iLi^i{®)/^' 
where Iat is the x 1 all one vector 

For the first objective function in ([8]l, as is decreasing 

and Schur-convex, the objective function satisfies ||9l 

^PM{'^N - e)) > (^1n - CEti^^i®)/^) ® 1^) ' 

=9[A(©)] 

(37) 



with equality holds if and only if [0];,^ = X]i=i^i(®)/-^- 
Therefore, must have the following structure ||9l 

= QFdiag(A(0))Qg. (38) 

where Qf is a unitary matrix such that has identical 
diagonal elements. 

Based on the definition that V'i(*) is a decreasing and 
Schur-convex function, based on A.6.Lemma and A.S.Lemma 
in ifTTl it can be directly proved that g(A(0)) is a decreasing 
and Schur-concave function of A(0). 

Obj 2: Notice that for the positive semi-definite matrix 
*MSE = Ijv - 0, d{lN - 0) -< A(lAr - 0) Furthermore 
ip2{') is Schur-concave, we have 

t/>2(cl(Ijv-0)) > V2(lw-A(0)). (39) 

=9lM@)] 

In order to make the equality in ( [39] ) hold, we need = 
Ai(0), which means that is a diagonal matrix. Therefore, 
we can write 

= lArdiag(A(0))lAr. (40) 

Since '02(*) is increasing and Schur-concave, based on 
A.6.Lemma and A.S.Lemma in 111] it is obvious that 
■02(lAf — A(0)) is decreasing and Schur-concave with re- 
spective to A(0). 
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